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In this paper, we consider Bianchi type III and Kantowski-Sachs space- 
times and discuss the behavior of time-varying constants G and A by us- 
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technique, we consider usual, modified, and completely modified matter 
collineation equations while studying the behavior of these constants. Fur- 
ther, we reduce the results for dust, radiation, and stiff fluids. We find that 
A is a decreasing time function while G is an increasing time function. This 
corresponds to the earlier results available in the literature for other space- 
times. Further, we find that the deceleration parameter attains a negative 
value, which shows that the expansion of the universe is accelerating. 
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I. INTRODUCTION 



The time- varying behavior of the gravitational and the cosmological con- 
stants, i.e., G and A, have been among the most controversial issues in cos- 
mology. Some authors reveal the idea that both G and A can be considered 
as non-constant, i.e., G = G{t) and A = A(i), coupling scalars, while solving 
the Einstein field equations (EFEs). The use of G and A in the EFEs plays 
a significant role in cosmological, astronomical, and quantum phenomena. 

Improved astronomical techniques indicate that a non-zero value of A 
is required to discuss distant supernova. Analysis of high red-shift super- 
nova shows that the universe may be accelerating [1-3] due to the presence 
of some type of vacuum energy. Some people relate this vacuum energy 
density to a non-zero cosmological constant. This indicates that the cos- 
mological constant is a cause of the expansion of the universe and plays an 
important role in the evolution of the universe. In gravitational collapse, the 
cosmological constant slows it down, which limits the size of a black hole [4]. 
However, there is a fundamental problem related with the value of A, which 
is assumed to be very small. This value is about 120 orders less than the 
magnitude of the vacuum energy-density calculated in quantum field theory. 
A phenomological solution to this problem is suggested by considering A as 
a function of time. In the early universe, A was large, but decreased with 
the expansion of the universe and the creation of photons [5] . 

The gravitational attraction, explained as a result of the curvature of 
spacetime, is proportional to G. According to quantum theory, the non- 
homogeneity of the gravitational field causes G to change rapidly in small 
intervals of time [6] . The EFEs with G and A are given by 

Rab-lR9ab-Ht)9ab = ^^-T ab , (a, b = 0, 1, 2, 3) , (1) 
— ■ c 

where R a b, g a b, and T a b represent the components of the Ricci, metric, and 
energy-momentum tensors respectively, and R is the Ricci scalar. 

The time variation of G was originally raised by Dirac [7]. He proposed 
that the gravitational constant varied with age of the universe. Modern 
theories, like string theory and Brans-Dicke (BD) theory, do not necessarily 
require such a variation, but provide a natural and self-consistent frame- 
work for this variation by assuming the existence of additional dimensions. 
The time variation of G in these multi-dimensional theories has recently 
been studied, and their consistency with variable observational data for dis- 
tant supernova has been analyzed [8]. The variation of the gravitational 
constant is found to make distant supernova appear brighter. The recent 
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results of Shapiro et al. [9], based on an analysis of radar echo time de- 
lays, have set an experimental upper limit on the possible time variation of 
the gravitational constant as \G/G\ < 3 x 10~ 10 yr~ 1 , where the dot repre- 
sents a derivative with respect to time. Theoretical calculations by Dicke 
for zero-pressure Friedmann-type cosmologies, according to the BD theory, 
yield \G/G\ « lO -11 ^ -1 for a flat spacetime and \G/G\ « 3 x lO -11 ^ -1 
for closed spacetimes [10,11]. Although the time variation in G is extremely 
small at the present epoch, Dicke has shown that there exist early epoch 
solutions for which the energy density associated with the time variation 
is much greater than the matter energy density. Recently, a constraint on 
the variation of G has been obtained by using the Wilkinson Microwave 
Anisotropy Probe (WMAP) and the big bang nucleosynthesis observations, 
which comes out to be -3 x 10 _13 yr _1 < (G/G) today < 4 x 10 -13 yr _1 [12]. 

Bekenstein [13] and Bertolami [14] introduced models in which both G 
and A are time dependent. Several authors [15-17] studied the variations 
of G and A in the framework of flat Friedmann-Robertson- Walker (FRW) 
symmetries. This work has been extended [18-21] to more complicated ge- 
ometries like the Bianchi type I model, which is the simplest generalization 
of the FRW flat model, by using a perfect fluid. The same model was 
considered with viscous fluids [22-24] to discuss the time variations of A 
and G. Kalligas et al. [25] discussed the behaviors of these varying con- 
stants by using Lie method. Darabi [26] found that the time variations of 
these constants lead the vacuum energy density to be time dependent as 
p v = A(t)/8irG(t). Therefore, for an early universe, where A is so large 
and G is so small compared with their current values, the vacuum energy 
is huge. At the present status of the universe, however, the vacuum en- 
ergy is vanishing due to the time variations of both A and G. Belinchon 
and Dvila [27] discussed time-varying constants in different spacetimes by 
using different symmetric techniques. In recent papers [28, 29], the same 
author analyzed the behaviors of time- varying G and A for a Bianchi type I 
model and made comparison of different techniques, including self-similarity, 
matter collineations, kinematic self-similar similarity, and the Lie method. 

In this paper, we extend Belinchon's work to Bianchi type III and Kantowski- 
Sachs spacetimes. We shall use two symmetric techniques, i.e., matter 
collineations (using energy-momentum tensor for a perfect fluid) and kine- 
matic self-similarity (first and second kinds). The scheme of this paper is 
as follows: In Section II, we shall write the field equations with relevant 
quantities for Bianchi type III and Kantowski-Sachs spacetimes. Section 
III is devoted to a study of the behaviors of the time-varying constants by 
using kinematic self-similarity technique. In Section IV, we use the matter 
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collineations technique to investigate the behaviors of G and A. The last 
section will provide a summary and a discussion of the results obtained. 



II. BIANCHI TYPE III AND KANTOWSKI-SACHS 
SPACETIMES 

Bianchi type III and Kantowski-Sachs spacetimes are spatially homogeneous 
spacetimes that admit an abelian group of isometries G3 acting on a spacelike 
hypersurface. These are generated by spacelike Killing vectors £1 = d r , £2 = 
8q and £3 = d$. In co-moving coordinates, the metric representing these 
spacetimes is written as [30] 

ds 2 = c 2 dt 2 - A 2 (t)dr 2 - B 2 (t){d6 2 + f 2 (9)d(t> 2 ), (2) 

where A and B are arbitrary functions of t while f{6) is defined as 

f{9) = sinh# corresponding to Bianchi type III spacetime, 
f{6) = sva.6 corresponding to Kantowski-Sachs spacetime. 

These metrics represent anisotropic generalizations of the open and closed 
FRW models, respectively. 

The energy-momentum tensor for a perfect fluid is given by 

Tab = {p + p)u a u b - pg a b, (3) 

where u a is the four- velocity and in co- moving coordinates, it is defined as 
follows: 

u a =Q, 0,0,0j withu a u a = l. (4) 

A perfect fluid can be characterized by a dimensionless number k given by 

P 

k = — or p = kp. (5) 
P 

This is called the equation of state and represents a dust fluid for k = 0, 
radiation for k = 3, and stiff matter for k = 1. Using Eq. ((H), the EFEs of 
Eq. ([T|) lead to the following three equations: 

2AB c 2 f" B 2 8vrG . , 

-AB-Wf + W = ^ + Ac ' ^ 

2B c 2 f" B 2 8TTG, 2 

+ -2 = 2-/cp + Ac 2 , (7) 



B B 2 f B 2 c 2 
A B AB 8vrG 
A + B + AB ~ 6 



+ o+TU = -kp + Ac 2 . (8) 
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Here, the dot denotes the time derivative, and the prime denotes the deriva- 
tive with respect to 6. 

The time derivatives of G and A can be related by the Bianchi identities 
as follows: 

(Rab ~ \Rgab)' b = {^T ab + Ag ab y b . (9) 

Simplification of this expression by fixing a = and varying b = 0,1,2,3 
yields 

p + p(l + k) \- A +2^- 



8vr ■ , 8ttG 
^p + A = — 



A B 



(10) 



The conservation law of the energy-momentum tensor of matter field, T% b = 
0, gives 

p + p{l + k)(^ + 2^j =0. (11) 
In view of this equation, Eq. (|10p implies that 



A = -^Gp. (12) 

It is mentioned here that for the possibilities a = 1,2,3 and b = 0,1,2,3, 
Eq. ([9]) is satisfied identically. 

We define Hubble's parameter as an average expansion of the universe 
as follows: 

H = j + 2^ = H 1 + 2H 2 , (13) 

where 

~ A' H2 ~B- 
Consequently, the deceleration parameter is defined as 

Using Eq. ^ in Eq. §\A§, it follows that 

p + p{l + k)H = 0. (15) 

In the next two sections, we use a Bianchi type III spacetime to discuss the 
behaviors of G and A. However, we include the discussion of the behaviors 
of G and A for Kantowski-Sachs spacetime in the last section. 
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III. TIME- VARYING BEHAVIOR OF G AND A 
BY USING A KINEMATIC SELF-SIMILARITY 

TECHNIQUE 



In general relativity, self-similarity can be defined by the existence of a 
nomothetic vector field. Cahill and Taub [31] were the pioneers to introduce 
the concept of self-similarity corresponding to homothety. Carter and Hen- 
riksen [32] introduced the concept of kinematic self-similarity as a natural 
generalization of the homothetic case. A kinematic self-similar vector field 
£ satisfies the following conditions [33]: 



where a and 5 are constants and h ab = — u a Ub is the projection tensor. 
Kinematic self-similarity (KSS) can be classified into first and second kinds 
by using a scale independent ratio j, referred to as the similarity index. 
The ratio j = 1 leads to self-similarity of the first kind, also known as ho- 
mothety, while the ratio j ^ 0, 1 indicates self-similarity of the second kind. 
Here, we discuss the behaviors of G and A by using these kinds. 

1. Kinematic Self-similaity of the First Kind 

For a = S = 1, the KSS vector field yields 



In this case, £ is known as a homothetic vector field. This gives the following 
system of ten equations: 



£^u a = au a , 
£^h a b = 25 h ab . 



(16) 
(17) 



9ab,d C + 9ad% + 9bc£% = ^9ab- 



(18) 



c 2 ^ - A% 
c 2 £° 2 - B% 
c 2 £°3-£ 2 sinh 2 ^ 3 
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(19) 
(20) 
(21) 
(22) 



i. 



(23) 

(24) 
(25) 



A 2 ^ + B 2 smh 2 e^ 



0, 
0, 



i 



(26) 
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B 2 ^ + B 2 sinh 2 9£% = 0, (27) 

(|)£°+£ 3 3 = 1. (28) 

Here, the derivatives with respect to 0, 1, 2, and 3 denote the partial deriva- 
tives with respect to t,r,6, and (ft, respectively. Solving Eqs. (fl~9j) . ([23]) . 
(|26p . and (|28p simultaneously, we obtain 

£° = t + f (r,e,<t>), (29) 
= (l-^°)r + /i(W), (30) 
e 2 = (l-§£> + / 2 (*,r,0), (31) 

e 3 = (l-§£> + /3(i,r,0), (32) 

where /i(t, 0, </>) — /s(i, r, 0) are functions of integration. Substituting these 
values of £ in Eqs. ([23]) and (|26|) . respectively, we obtain 

2AA/ o , 1 (r,e,0)r = O, (33) 
2BBf 0i2 (r, 9,^)6 = 0. (34) 

Since A,B,r, and # 7^ 0, the above equations yield the following four cases: 

(0 A = = B, (ii) A = 0, 5^0, (m) i^0, S = 0, i, B + 0. 

The cases (i)-(iii) do not provide the behaviors of G and A as they vanish. 
In the following, we discuss the case when A, 5^0. 
In this case, Eqs. (f33j) and (f34"j) yield 



/o(r,M) = / o (0), 
which, in view of Eq. (|28p . finally gives 

/o(0) = c , 

where Co is an arbitrary constant. Thus, £° takes the form 

f = t + c . (35) 
When we make use of this value of £° in Eqs. (|30p - (|32p . it follows that 

f 1 = (l-j(t + co))r + fi(t,9,<l>), (36) 
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f = (l-§(i + eo))0 + / 2 (i,r,£), (37) 

£ 3 = (l-^(t + c ))<t> + Mt,r,9). (38) 

Now replacing these components of £ in the above system of equations, 
the functions of integration fi(t, 0,(f>) — f3(t,r,9) reduce to arbitrary con- 
stants ci, c 2 , and C3 respectively along with the following ODEs as necessary 
and sufficient conditions: 

AA + (t + c )(AA - A 2 ) = 0, (39) 
BB + (i + cq)(BB — B 2 ) = 0. (40) 

Thus, the homothetic vector field becomes 
C = (t + c )d t + [(l-(t + c )j)r + c 1 }d r + [(l-(t + c )^)e + c 2 }d e 

+ [(l-(t + co)|)0 + c 3 ]^. (41) 

Using Eq. ([n?]) , the solution of the ODEs in Eqs. (]3"9"]) and (|4T)]) . respectively, 
yield 

A = A (t + c ) Ql , B = B (t + c ) Q2 , (42) 

where Aq, Bq, a%, and a 2 are positive constants (for physical reasons) of 
integration. Making use of these values of A and B in Eq. (|13p . the Hubble 
and the deceleration parameters turn out to be 

H = (a 1 + 2a 2 )(t + c )-\ (43) 

Q = * -1- (44) 
«i + 2a>2 

Consequently, Eq. (fT5|) yields 

P = Po(t + c )- (1+fc)(Ql+2a2) , (45) 

where is a constant of integration. 

Solving Eqs. © and ((7j) simultaneously and then using Eqs. (132]) and 
150. we have 

G = g + fc) («1"2 ~ «2(« 2 - l))(t + Co )-2+(tti+2a 2 )(l+fc) ) (46) 

1 2 
A = -^[2aia 2 + a 2 2 - - ; , (aia 2 - a 2 2 + a 2 )](t + c )~" 



c 2L 1 + Jfc 

1 



4(* + co)^ 2 . (47) 



'0 



S 



Using the values of G, and A and Eq. ([35]) in Eq. ©, it follows that 

c 2 

a 2 (a 2 - 1) + ai«2 + ai(«i - 1) = 2a 2 (a 2 - 1) + a 2 2 - -=r^(t + c )~ 2a2+2 . 

Bo 

This shows that ct\ will be constant only if we choose a 2 = 1. Thus, we 
have 



For cti to be real, it is necessary that 



Bo 2 ' 



B 2 >c 2 => B >c or B < -c. 
In other words, we can say that 

B G (-oo, -c] U [c, do). (49) 
For Bq, to be positive, our interval of interest is only where Bo £ [c, oo). 

i. Behavior of G 

Using ot 2 = 1 in Eq. (JM|), G becomes 

G = G (t + c r +(ai+2)fc , (50) 

where 

Go = 2CV M . (51) 
We note that Eqs. (05} and dSDJ yield 

Gp^{t + co)- 2 . (52) 



For Go > 0, i.e., a\ = yl — -^y, the behavior of G can be discussed as 
follows: 

„. . . -2k 

G is increasing 44> ai > 



(1 + fc) 



c 2 -2fc , . 

* V 1 <53) 

When we take square of both sides, the above inequality yields the following 
two cases: 



(i) Here, the inequality in Eq. ([53]) implies that 



c 2 4k 2 „ / c 2 2fc 

1 - — ^7 > if W 1 - r=-n > 



Thus, 



5 2 (1 + A;) 2 V V (1 + fc)' 



n ■ ■ d 2 / c(l + fe) 

G is increasing 44> iso > 



where 



V1 + 2&-3A; 2 , 
-B G (-oo, -6) U (b, oo), (54) 

C(1 + &) V *€(-!,!); (55) 



Vl + 2A; - 3/fc 2 ' v 3' 

6 becomes imaginary or infinite for all other values of k. One can easily 
verify that b > c V fe € (— |, 1). 

(ii) The inequality in Eq. (f53l) implies that 



c 2 4A; 2 „ / c 2 2ifc 

1 - — n < if \ / 1 - tn? < 



Thus, 



5o 2 (l + *0 2 V #o 2 (1 + fc)' 



^ • • d 2 / C(l + fe) 

G is increasing 44> Bq < 1 



y/l + 2k-3k 2 , 
B G (-6,-c)U(c,fe), (56) 



where 6 is the same as defined in Eq. (|55p . Hence, from both the cases, 
we can conclude that G increases for all -Bo G ^ + \(0,c) while G becomes 
constant at b. It is mentioned here that G always vanishes at Bq = c. 

Now, we discuss the behaviors of G in the dust, radiation, and stiff fluid 
cases. For dust, we take k = 0, and the behavior of G is the following: 

G is increasing 44> -Bo G 5ft + \(0, c), 
G vanishes 44> Bq = c. 

For the radiation case, we have k = 1/3; thus, 

2c 

G is increasing 44> -Bo G 3? + \(0, —7=), 

v3 

2c 

G is constant 44> -Bo = — y=. 

v3 

For the stiff fluid, G increases for all values of -Bq G 5R + \(0, c). 
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ii. Behavior of A 

To discuss the behavior of A, we substitute «2 = 1 in Eq. (|47p so that A 
becomes 



A 



(2ai + !)(! + k)B 2 -c 2 (l + k)- 2 ai B 2 



c 2 B 2 (l + k) 



(t + c )- 2 = A (t + c )- 2 , 
(57) 



where 



if 


Ao 


< 


and 


t > c , 


or 


Ao 


> 


and 


t < c , 


if 


Ao 


> 


and 


t > c , 


or 


Ao 


< 


and 


t < c , 


if 


Ao 


= 


or t - 


-> oo. 



. (2a 1 + l)(l + /c) J B 2 -c 2 (l + k)-2 ai B 2 
A ° " c 2 B 2 (l + k) (58) 

is a constant. From this equation, we can discuss the behavior of A as 
follows: 

A is increasing 

A is decreasing 

A vanishes if Ao = or t — > oo. (59) 

It is noticed that Ao > if 

(1 + k){2 ai B 2 + (B 2 - c 2 )) > 2 ai B 2 . 

Using Eq. (08]), this inequality becomes 

(1 + fe)(2ai + a\) > 2ai, (60) 
-2k 

* ai> TTI' 

which is the same condition as given by the inequality in Eq. (|53|) . Hence, 
we finally obtain 

A >0 V B g K + \(0,c). (61) 

Similarly, Ao < if 

(1 + k)(2a x + q 2 ) < 2ax, 

which yields a contradiction. Further, Ao vanishes at Bq = b. 

For different types of fluids, we can discuss the above conditions on A 
as follows: In the dust and stiff fluid case, Ao is always positive for any 
value of Bq E 9? + \(0, c) while in the radiation case Ao is positive for all 

s e^ + \(o,^§). 



n 



2. Kinematic Self-similarity of the Second Kind 

For a = 5 ^ 0, 1, the definition of KSS vector field yields 



c 2 ^ - A 2 ^ 

9 v-D ^9 

c 2 e° 3 - # 2 sinh 2 #£ 3 


= 2ac 2 , 

= o, 
= o, 
= o, 


(62) 
(63) 
(64) 
(65) 






(66) 


A 2 ^ 3 + B 2 sin 2 0£ 3 i 


= o, 
= o, 


(67) 
(68) 




= 26, 


(69) 


e 2 3 + sinh 2 0$ 


= o, 


(70) 


(f)*° + $ 


= 25. 


(71) 



Solving this system of equations simultaneously by adopting the same pro- 
cedure as in the first kind of KSS, it follows that 

£ = (at + p)d t + [{6-j(at + p))r + c 5 }d r 

+ [{6-^(at + (3))e + c 6 ]d e + [(5--{at + (3))<f> + c 7 ]dt. (72) 

D A 

A = A (t+P) a \ B = B (*+-) aa , (73) 
a a 

H = (a 1 + 2a 2 )(t + -)- 1 , (74) 

a 

p = P0 (t + P)-v+m«i+*») j (75) 

a 

where Aq, Bq, ati, a 2 , and po are constants of integration. 
Similarly, G and A can be found to be 

G = G (t + -r i+{Ql+2)fc , (76) 
a 

A = A (t+^)" 2 , (77) 
a 

where Go and Ao are the same as given in the first kind. Obviously, Eqs. 
(|76"j) and ([77]) reduce exactly to Eqs. (foTJj) and (|57p for cq = -; hence, the 
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behavior of G and A will be the same as in case of KSS of the first kind just 
by replacing cq = @- in the corresponding constraint equations. 

IV. TIME- VARYING BEHAVIOR OF G and A 
BY USING THE MATTER COLLINEATION 

TECHNIQUE 

The energy-momentum tensor T ab represents the matter part of the field 
equations in Eq. ([TJ) . This enables us to understand the physical struc- 
ture of spacetime. Symmetries of the energy-momentum tensor provide the 
conservation laws on matter fields. Matter collineation is defined as 

£(T ah = 0. (78) 

This does not give any information about the behaviors of G and A. How- 
ever, it can be modified to obtain the behavior of G as follows [28]: 

M 8 -^T ab ) = 0. (79) 

If we introduce A, we get a complete modification of matter collineation, 
Eq. ([TED [28]: 

£ ^^m Tab + k{j . )gah)=Q {m) 
This equation helps us to discuss the behaviors of time-dependent G and A. 
1. Matter Collineations 

Using the components of a perfect fluid energy-momentum tensor, the sys- 
tem of matter collineation (MC) equations yields 

(81) 
(82) 
(83) 
(84) 
(85) 
(86) 
(87) 



(pc 2 ),oe + 2pc 2 e°o = 


= o, 


pc 2 $+pA% = 


= o, 


pc 2 C° 2 + P B 2 ^ = 


= o, 


pc 2 i% + P B 2 sinh 2 #£ 3 = 


= o, 


(pA 2 ) fi e + 2pA% = 


= o, 


pA 2 i] 2 +pB 2 i 2 l -- 


= o, 


pA 2 ^+pB 2 smh 2 e^ x = 


= o, 


(pB 2 ) fi e + 2pB 2 £ 2 2 = 


= o, 
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pB 2 £ 2 3 +pB 2 swh 2 6£,% = 0, (89) 
{pB 2 sinh 2 6) fi f + 2pB 2 sinh 2 6£% = 0. (90) 

It follows from Eq. dgTJ) that 

£° = /8(r,M>- 1/2 , (91) 
where fs(r,9,(p) is a function of integration. Similarly, Eq. (j85|) yields 

(^ + 24K + 2^ = 0. 
p A 

Using Eqs. ([5]) and (f9Tj) . this equation gives 

f 1= ^(^-l) r + M*,e^), (92) 

where fg(t,9,4>) is an integration function. Similarly, by solving Eqs. (|88p 
and (|9U|) and then integrating w.r.t. and </>, respectively, we obtain 

f = e°^-|^ + /io(t,r,0), (93) 

e = £°(|S-§^ + /ll(i,r,0) 5 (94) 

where fio(t,r,(f>) and fu(t,r,0) are functions of integration. 

When we solve the above system of ten equations simultaneously by using 
these values of £, as done in the KSS technique, the functions of integration, 
fs{r,6,4>) — fu(t,r,8), reduce to arbitrary constants, which are termed as 
C8,cg,cio, and en, respectively. Thus, the final form of the vector field, in 
terms of £° = csp~ 1 ^ 2 , can be written as 



c = fd t + Ke° - jfy + C9 r + [(e° - §e> + Cl0 ]<% + [(e° - + Cll ] 

along with the following constraint equations 



(95) 



These two equations reveal that either p = 0, i.e., a dust case, or p ^ 0, 
i.e., a perfect fluid case. The dust case is not interesting as the time- varying 
behaviors, of G and A cannot be discussed for this case by using MCs. 
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In the perfect fluid case, Eq. ([96]) implies that 

A = —— e ai I VP dt + a 2 ^ £> = _L e bl I Vpdt+ki (97) 

where a\,ai,bx, and 62 are arbitrary constants of integration. For these 
values of the metric functions, the EFEs give 

G ~ 8^(1 + k) [ 7 ~ 7 Cl3 7^ ~W 1 m 

A = i[(6 1 //2_A)2_ ( ^ e -6 1 /^ + 6 2) 2 
c z 2p p L ' z 

2k pp^ 2 p 2 p- 2 

+ TTk {b2p - Cl5 — + — ) 

+ —k^^J^-"^-^^ (99) 

These equations show that the time-varying behaviors of G and A cannot 
be discussed unless p is given. We assume the following two cases to discuss 
the behaviors of G and A: 

(i) p(t) = po (constant); 

(") />(*) = pTio 7 ' where ° = i and b = %■ 

Case (i) When we substitute p(t) = po (constant) in Eqs. (|98|) and ([99]) . it 
follows that 



and 



where 



A= — e - 2bl(t+b ^ +A , (101) 
Po 



a 1 r/i 1/2n2 2A; 2 61 

Ao = ^[(M ) +rri^ + TTI (ci3P0 + ^ )] - 

We see from Eqs. (jlOOp and (|10ip that G always remains constant while the 
behavior of A is as follows: 

A is increasing 44> b\ > 0, 
A is decreasing 44> 61 < 0, 
A is constant 44> b\ = 0. (102) 
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Case (ii) In this case, we get £° = t + cq. Using this in Eq. (|95p . we obtain 
the same homothetic vector field as in case of KSS of the first kind given in 
Eq. (|4ip . It is mentioned here that this homothetic vector field satisfies the 
relation 

£ HO (T a b) = (103) 



with the same ODEs as given by Eqs. (|39p and (|57p . Thus, the behaviors 
of G and A will be the same as in the case of KSS of the first kind. 

2. Modified Matter Collineations 

The modified MC equations, Eq. (|79f) . can be written as 

i°GT ah + G(T abjC e + Tactf, + TU%) = 0. (104) 
The corresponding system of MC equations yields 

= cMGp)' 1 / 2 , (105) 

e = Z°($-f)r + ci8, (106) 

e = e°(|-|)^c 19) (107) 

e = *° (108) 

Here, Ci7,ci8,cig, and C20 are arbitrary constants of integration satisfying 
the following ODEs 

pdt[e{1 2 { -p + ^ ) + J )] = °' (109) 

^ [e ° ( i ( p + § ) + § )] = °" (110) 
For p ^ 0, Eqs. (fTMl) and (fTTUj) provide 

= - c c 21 J dt+C 2 2 g _ _ £ C23 / V 7 ^ dt+C 2 4 (Hi) 

For these values of A and B, the EFEs yield 

- 
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- {c2,{Gp)^ 2 - \{ P - + ^)} 2 - \c 2h (Gp + Gp){GpT 112 

A = ^[{c 2 5(Gp)^-I(^ + |) } 2_ { -^ e --/^ +C25}2 

+ ^Z{^(GP) 1/2 ~ \£ + %)H<*(Gp)W - \{t + £)} 
1 + « 2 p G 2 p G 

1 + « 2 p G 2 

- -(^- p - + ^--)M (113) 

It is obvious from Eqs. (|112p and (|113p that the expressions of G and A are 
too complicated to discuss. Thus, we assume £° = t + Co, i.e., G(t)p(t) = 
( at + b yi (as found in homothetic case), such that the vector field becomes 
homothetic and satisfies the equation 

£ho M TaS ) = (114) 



along with the same ODEs as given by Eqs. (|39p and (|40p . This corresponds 
to the KSS of the first kind. It is worth mentioning here that this case leads 
to 

% + - = tttN °P ~ (* + c o)~ 2 - ( 115 ) 
G p (i + c ) 

which is exactly the same relation for the product Gp as found by using the 
KSS technique. 



3. Completely Modified Matter Collineations 

This is the case in which the variation of A is included, as well. Here, we 
consider Eq. (|80|) and apply the definition of the Lie derivative to obtain 

eGT ab + G{T ab , c C + T a £% + TU%)] 

+ [M°9ab + H9ab, c e + dac^b + 9bcC%) = 0. (116) 

The system of equations here leads to the following solution: 
£0 _i 



{Gp + pG) + A 
H (Gp)+A 
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giving 



£° = c 26 (r,0,0) 
The remaining components of £ are 



+ A 



-1/2 



e 3 



(Gp + pG) - A 2A 
(Gp + pG) - A 25 

+ IF 



fr(Gp) - A 



-1 



^(Gp + pG) -A 2B_ 
H{Gp)-A + B 



Z°8 + c 28 (t,r, <j>), 
£V + c 29 (i,r,0). 



(117) 

(118) 
(119) 
(120) 



The following ODEs become the necessary and sufficient conditions: 



t 
d t 



^(Gp + pG)-A i 2A\{8irG P \~ 1 / 2 



8tt 
1* 



(Gp)-A + A ) \ c 4 
(Gp + pG)-A 2B\ (8ttGp 



+ A 



8tt 



(Gp)-A 



+ 



+ A 



-1/2 



With the equation of state and then replacing 



provide 



SwGp 



= 0, (121) 
= 0. (122) 
the above ODEs 



A 
B 



A ($ - A)" 1 ' 2 + e as I Vf+^ d \ 
B ($ - A)" 1 / 2 + e a4 / * 



(123) 
(124) 



where Aq, Bq, a^, and 04 are constants of integration. If we substitute these 
values of the metric functions in the EFEs, we get much more complicated 
expressions for G and A. However, if we assume £° = t + Co, then the 
following relation is satisfied: 



£ho 



/87rG(i) T ^ + 



0, 



yielding the metric functions in the form 



A = 4)(* + co) Ql (*-A)- 1/2 , 
B = B (t + c ) a2 ($-A)- 1 / 2 . 



(125) 



(126) 
(127) 



Here, ol\ and a 2 are constants of integration. When we use these values 
of the metric functions in the EFEs, we are again unable to discuss the 
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behaviors of the G and A due to their complicated expressions. We note 
that if we replace A by $ — 1, i.e., A = 87T ^ P — 1, then it also corresponds to 
the homothetic case of the KSS. 

It is interesting to note that if we write Eq. (|80p as 

£ i (^r^ T «*) = = £ d A (t)9ab), (128) 

which is a special case, we may get some insight. In this way, we reach again 
homothetic cases and obtain the following results: 

Gp^(t + c )- 2 and A = A (i + c )~ 2 , (129) 

where Ao is an arbitrary constant of integration. Further discussion on the 
behavior of A is the same as given in Eq. (|59p . 

V. SUMMARY AND DISCUSSION 

We have studied the perfect fluid Bianchi type III and Kantowski-Sachs 
spacetimes with time- varying constants G and A. Due to the time- varying 
nature of these constants, Bianchi identities, along with the energy-momentum 
conservation law, = 0, yield a time-dependent expression of the energy 
density. This expression helps us to define Hubble parameter and the decel- 
eration parameter. 

In the KSS, we studied the behaviors of G and A for the first and the 
second kinds. When we solve the system of ten self-similar equations simul- 
taneously, there arise two ODEs as necessary and sufficient conditions. The 
solution of these ODEs yields the metric functions A and B, which make 
the variational behaviors of G and A possible. We also discuss the dust, 
radiation, and stiff fluid cases. Further, we discuss the behaviors of G and 
A by using MCs only the case that corresponds to the KSS. 

In the KSS of the first kind, the metric functions take the form 

A = A (t + co) ai and B = B (t + c ) Q2 , 

where a\,a2,Ao, and B$ are arbitrary positive constants. It is worth men- 
tioning here that the physical situation is only possible if we assume «2 = 1- 
Then, G takes the form 

G = G (t + c ) ai+{ai+2)k , 
where Go is a constant given by Eq. (1511) . 
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For a Bianchi Type III metric, ct\ = J 1 — and the behavior of G 
depends on b, Eq. ([55]) . For the Kantowski-Sachs metric, ot\ turns out to 
be yjl + and the behavior of G depends on d = r/gpjjjj 2fc V G 
(-oo,-±)U(l,oo). 

Table 1. Behavior of G for both spacetimes. 



Spacetime 


G is increasing 


Constant 


Bianchi type III 


B e R+\(0,c) 


B = b 


Kantowski-Sachs 


Any value of £? G K + 


B = d 



The dust, radiation and stiff fluid cases for Bianchi Type III metric are given 
in Table 2. 

Table 2. Behavior of G for different cases of fluids. 



Behavior of G 


Dust case 


Radiation case 


Stiff matter 


G is increasing 


£ G»+\(0,c) 




B G K+\(0,c) 


G is Constant 


5o = c 







In the case of the Kantowski-Sachs metric, the dust and radiation cases are 
not physical while for stiff matter, G is increasing V Bq G In the KSS 
of the first kind, A is given by the following equation: 



A = A (t + c ) 2 , 

which yields the behavior of A given below in Table 3. In Bianchi type 
III and Kantowski-Sachs metrics, the possible values of Ao are discussed in 
Table 4. For various fluids, Ao reduces for Bianchi type III, as given in Table 
5. In Kantowski-Sachs metric, the solutions corresponding to the dust and 
the radiation cases are not physical while for stiff matter Aq > V Bq G 5i + . 



Table 3. Behavior of A. 



Behavior of A 


Value of Ao and t 


A is increasing 


Ao < and t > c$ or Ao > and t < Co 


A is decreasing 


Ao > and t > c$ or Ao < and t < Co 


A vanishes 


Aq = or t — > oo 



Table 4. Possible values of Aq. 



Spacetime 


A > 


A = 


Bianchi type III , 


B G »+\(0,c) 


Bo = b 


Kantowski-Sachs 


V B G 


B = d 



20 



Table 5. Possible values of Aq for different types of fluids. 



Cases 


Dust Fluid 


Radiation Case 


Stiff Matter 


A >0 


S G5R+\(0,c) 


5oer\(0,f 


B e5R+\(0,c) 


A = 


B = c 


R 2c 





We note that for the KSS of the second kind, the results coincide with those 
of the first kind as given in the above tables except that in the expressions 
for the metric functions, cq is replaced with a fraction of two constants 

Using the MCs technique, the behaviors of G and A are not straight- 
forward as in the case of the KSS due to the complicated metric functions. 
However, we have managed to discuss two particular cases depending upon 
p. Firstly, for p = po (a constant), G becomes constant while A varies 
as given in Eq. f)102|) . Secondly, the case for p = / at + b yi , corresponds to 
the homothetic case. In modified MCs, the behaviors of G and A could 
not be discussed generally, but we obtained a homothetic case by assuming 
~ (at+b)' 2 " Further, we obtained a relationship Gp « (t + cq) 2 , 
which was the same as obtained in the KSS technique. Similarly, in the 
completely modified MCs case, we again obtained a homothetic case by as- 
suming £° = (t + Co) and A = 8n ^[ p — 1. Further, when we re-interpreted 
the completely modified MC equations, given in Eq. ([80]) . by Eq. (|128p . we 
directly obtained a homothetic case giving A = Ao(i + Co) -2 . Consequently, 
the behavior of A turns out to be the same as given in Eq. (|59p with the 
KSS. It is mentioned here that the metric functions and vector fields are the 
same for both spacetimes. However, we obtain different behaviors of G and 
A due to slight changes in the EFEs. 

We found that the cosmological constant A turned out to be a time 
decreasing function for Ao > while the gravitational constant G was a 
time increasing function when Go > for all values of k > 0. It is worth 
mentioning here that our results verify the results obtained by Belinchon 
[29]. For these behaviors of G and A, the time-dependent vacuum energy 
density relation is also satisfied, according to which both these constants 
are changing in a reciprocal way [26]. Further, we found that for 02 = 1, 
the deceleration parameter attained a negative value, which showed that the 
expansion of the universe was accelerating. Thus, we can say that with the 
expansion of the universe, A is going to reduce [5]. 

We would like to mention here that the above mentioned time-varying 
behaviors of G and A can only be discussed in the homothetic case, i.e, 
the KSS of the first kind for both Bianchi type III and Kantowski-Sachs 
spacetimes. Moreover, we found that the vector field satisfying equation 
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£^9ab = 2gab also satisfied equation £(T a i, = 0. By modifying the MC 
equations in an appropriate way, we were able to find the same relationships 
as in the case of the KSS solution. 
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